JLAB-THY-02-59 



Radiative corrections and parity- violating electron-nucleon 

scattering 



S. Barkanova^, A. Aleksejevs^, and P.G. Blunden^'^ 

^Department of Physics and Astronomy, University of Manitoba, Winnipeg, MB, Canada 

R3T 2N2 

"^Jefferson Laboratory, 12000 Jefferson Avenue, Newport News, VA 23606 



Abstract 



Radiative corrections to the parity-violating asymmetry measured in elas- 
tic electron-proton scattering are analyzed in the framework of the Standard 
Model. We include the complete set of one-loop contributions to one quark 
current amplitudes. The contribution of soft photon emission to the asym- 
metry is also calculated, giving final results free of infrared divergences. The 
one quark radiative corrections, when combined with previous work on many 
quark effects and recent SAMPLE experimental data, are used to place some 
new constraints on electroweak form factors of the nucleon. 
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I. INTRODUCTION 



One of the areas of active theoretical and experimental interest is in understanding the 
strange quark contribution to nucleon electroweak form factors. The s quarks are a key 
to the overall sea properties, and their distribution is of particular interest for developing 
our understanding of low energy nucleon structure. At present, theoretical estimates of 
strangeness electric and magnetic form factors span a wide spectrum in both magnitude and 
sign [1]. 

The electromagnetic current is a pure SU(3)/ octet. By exploiting flavor symmetries one 
can separate two flavor structures: the isovector, (u — d), and the hypercharge, {u + d — 2s). 
There is no singlet component {u + d + s), so the electromagnetic current can not provide 
sufficient information to separate u, d, and s contributions. 

The offers a new flavor coupling to the nucleon proportional to weak isospin, which 
samples {u — d — s) in the light quark sector. Knowing matrix elements of {u — d), the Z*^ 
can be used as a probe to flnd strange quark matrix elements in the nucleon [2]. One can 
measure such direct characteristics of strangeness in the nucleon as the strange magnetic 
moment fig = x as, the strangeness charge radius = s'^r^s, and the strangeness 

analog of the axial charge 57755,, which can serve as an independent conflrmation of the 
quark spin fraction measurements. 

Electroweak properties of the nucleon can be studied by parity- violating electron scatter- 
ing at low to medium energies [3]. There the asymmetry factor coming from the difference 
between cross sections of left- and right-handed electrons can be measured. However, ex- 
tracting information of interest from the experimental data requires evaluating radiative 
corrections to electroweak scattering at the few percent level. 

Electroweak radiative corrections to intermediate energy, parity non-conserving semi- 
leptonic neutral current interactions have been addressed previously [4-8]. In these works, 
radiative corrections are constructed from the underlying fundamental weak interaction be- 
tween electron and quarks. Broadly, such corrections are denoted as being either one-quark 
or many-quark effects. The one-quark corrections involve the interaction of the electron 
with a single quark. The many-quark contributions involve two or more quarks, and include 
effects due to an intrinsic weak interaction in the nucleon (e.g. the anapole moment). In 
this work, we restrict our considerations to the one-quark contributions. 

The Tisc of computer packages FeynArts, FormCalc, and LoopTools [9] is of great 
assistance in facilitating our calculations, allowing us to include the full set of one-loop con- 
tributions (there are several hundred Feynman diagrams, which are laborious to calculate by 
hand), and to retain the momentum-dependence of the amplitudes, for example. Our treat- 
ment of hadronic model dependencies (e.g. kinematics) is slightly different than previous 
work. In addition, we have treated infrared (IR) divergences in the one-loop amplitudes by 
including bremsstrahlung contributions for soft photon emission. While such a prescription 
is not completely satisfactory, and would be better handled by also accounting for hard 
photons and the particular detector setup of a given experiment, it nevertheless sets the 
scale for uncertainties of this origin. 

The article is described as follows. In Sec. II we give the formalism for parity-violating 
electron scattering, including deflnitions of the relevant currents, couplings, form factors, 
and kinematics. Section III outlines the one-loop calculations, including the prescription for 
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including soft photon emission. Results for the radiative corrections are given in Sec. IV. 
Finally, our results are combined with other calculations of many-quark effects to discuss 
the implications for the SAMPLE experiment. 



II. FORMALISM 

A. Parity violating electron scattering 

At tree level, the electron-nucleon scattering amplitude consists of two terms, Ai'^ 
(Fig. la) and Ai^ (Fig. lb). The first term is parity-conserving 7-exchange, and the sec- 




FIG. 1. Tree level electron-nucleon scattering. 



ond has a parity- violating contribution from Z°-exchange. The differential cross section is 
proportional to 



da 



= |M^|' + 2Re{(A^^)*M^} + |M^ 



(2.1) 



The purely weak term 



is very small compared to the other terms and can safely 

be neglected. The electromagnetic-weak interference term 2Re |(7Vf^)* A/f^j contains the 
physics of interest. This term can be extracted from the parity- violating asymmetry, formed 
as the ratio of helicity-dependent to helicity-independent cross sections: 



dan ~ da I 
daji + dai 



(2.2) 



Here, a^ and a^ are cross sections for right- and left-handed electrons, respectively. This 
ratio, measured in a number of experiments, is of the order 10~^ at low energies. 
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For elastic electron- nucleon scattering, the asymmetry ratio is given by [7,10] 



A 



4\/27rQ; 



sGlGl + tG 



7 rz 



4 sin' ew)e'GlGl 



X 



N 



e{Giy + T{Gir 



(2.3) 



where Q' > is the four-momentum transfer. and G]^ are the electric and magnetic 
vector form factors of the nucleon associated with 7-exchange, Gf and Gf^ are the similar 
parameters for Z°-exchange, and is the axial vector form factor. Kinematic parameters 
T, e and e' are defined as 



4M2,' 



e — 



2 ' 



1 + 2(1 + r) tan 



(2.4) 



The electromagnetic vector form factors of the proton G^'^ and Gjj[ arc well measured, 
and some constraints can be put on G^" and G]^" for the neutron [11]. The expressions for 
the weak form factors of the nucleon, Gf and Gfj, and the axial vector coupling G^, are 
given below. The factor (1 — Asin^ 9w) — 0.1 suppresses the term containing G^, and makes 
the contribution of higher order processes more apparent. 

The neutral weak vector form factors for the proton and neutron can be expressed in 
terms of the electromagnetic form factors of the proton and neutron, plus a contribution 
from strange quarks [12]. We choose to write this in the form 



:i-2sin' 9 



w 



1 + R 



V 



-2 sin' 9 



w 



1 + 



r<T=o 



1 + 4°^ 



EM- 



(2.5) 



Here 7?^=°, 



pT=l 



and Ry^ are the isoscalar, isovector, and isosinglet weak radiative correc- 
tions describing the contribution from weak vector couplings beyond tree level, respectively. 
Strong isospin ra is equal to -|-1(— 1) for the proton (neutron). The isoscalar and isovector 
electromagnetic form factors are taken to be the linear combinations 



^E,M — ^E,M "T ^E,M^ 

gT=^ 



'E,M 



^EM ~ ^ 



EMi 



so that at Q' = 0, Gi=°(0) = G\=^{^) = 1, and G^=°(0) = 0.8797, G^=i(0) 
The strange quark form factors arc undetermined, but take on the values G|;(0) 



(2.6) 

4.709. 

and 



^m(O) — l^s: where fig is the contribution of strange quarks to the magnetic moment of the 
nucleon. 

The axial vector coupling is conveniently written in terms of octet matrix elements [7] 

as 
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1 + 



+ 



1 + i??) 



A- 



(2.7) 



The isovcctor axial form factor ^ is determined from neutron beta decay as ^(0) = qa, 
with = 1.2670 ±0.0035. The second term involving the 5*^7(3) isoscalar octet form factor 
is not present at tree level. Prom a least squares fit to hyperon beta decay, we have 
determined = 0.169 ± 0.009, which is consistent with results of polarized deep inelastic 
lepton scattering [13]. In our calculations (see section IV) we find < l-R^^"^!, so we 

expect overall that the second term in Eq. (2.7) is suppressed relative to the first term. The 
axial strange form factor G\ is extracted from polarized deep inelastic lepton scattering 
as G%{Q^ = 0) = -0.086 ± 0.024 [13]. Due to the unknown dependence we take 
G^(Q2 = 0.1) = -0.086 ± 0.086 in our analysis. 

Additional information on the axial vector form factor can be provided by parity- violating 
quasielastic scattering from deuterium. In the simplest impulse approximation, the asym- 
metry ratio for the deuteron can be written as the incoherent sum of neutron and proton 
contributions 



Ad 



-GfQ' 
4^2- 



TTQ! 



(2.8) 



where Np{Nn) is the numerator and Dp{Dn) is the denominator in Eq. (2.3) for the proton 
(neutron). 



B. Radiative effects 

As was indicated in the previous section, extracting form factors Gf and G^ is based 
not only on experimental data, but also on theoretical calculations of radiative corrections 
for parity-violating scattering. We calculate the radiative corrections to tree level for the 
electroweak interaction in the scattering processes e + p — > e + p and e + n^e + n hy 
modeling the nucleon as a collection of quasi-free quarks carrying some fraction x of the 
nucleon four-momentum. 

If we assume that during the scattering process the electron interacts with only one 
quark, we can split the problem of electron-nucleon scattering into calculations of Feynman 
graphs for processes like e{pi) + q{p2) ^(ps) -|- q{p4), with q representing {u, d, s) quarks. 
Starting from the fundamental coupling of an elementary fermion to the photon or to the 
Z^, one can arrive at the general form of electromagnetic and weak invariant amplitudes 
[7]: 

where /'^ (l^^) and ( J^^) are leptonic and hadronic vector (axial vector) currents, respec- 
tively, and Qf is the electromagnetic charge number of the fermion. 
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The weak coupling has been expressed in terms of the set {a, Mw, Mz) rather than 
Gp, noting that sin^ Ow = 1 - M^/M| in our renormahzation scheme. The vector and 
axial-vector "charges" of the fermion, gy and g-^, are defined as 

g{, = 2Ti - AQfSin^ (2.10a) 
g{ = -2r/, (2.10b) 

with Ti = +\{-\) for the upper (lower) member of the fermion doublet. 
Explicitly, we have for the electron and the {u, d, s) quarks: 

g^ = -l + Asiii^ew, 9%^+^, 

g-=+l-^sm'9w, 9l = -l, (2.11) 
4 

a.s 1 I • 2 /I d.s I 1 

9v =-l + gSm Ow, 9a = +1- 

For the case of electron scattering, leptonic vector and axial-vector currents are Dirac cur- 
rents with electron spinor Ug. 

It' = u^Yue, (2.12a) 

F^ = Uel^^^Ue. (2.12b) 

Hadronic currents J^^^ are hadronic matrix elements of the electro-magnetic, vector, and 
axial-vector quark current operators: 

J^^{n\j,\n), (2.13) 
(J, = j;, J'^ or J% ■ \N) = \p) or |n)), 

1 

^H9v"'^qllxUq, (2.14) 
9 

•^5 = ll9A^ql^lbUq. 

q 

Nominally, summation must be done over all quark flavors, but it is sufficient at the 
momentum scale of interest to include only light quarks, -u, d, and s. In order to be parity- 
violating, exchange must involve either V{e) x A[q) or A[e) x V{q). Thus, it is convenient 
to express the amplitude for the parity-violating (PV) part of the electron-quark {eq) scat- 
tering in the general form 



■2 



+S^«'Ue(P3)7''75Me(Pl) " Mg(P4)7M«<z(P2), (2.15) 
= Aeq jeq , peg jeq 

— ^ 'JvA + -'Ay- 
Here, and B'^'^ are functions of the Standard Model parameters, n— point tensor coeffi- 
cients, and kinematics. From this point we define one-loop radiative correction to tree level 
electron-quark scattering as follows: 
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p9 — rad _ -^rad /r) 

-"■A — Aeq 1 — r>eq i K"^-^^) 

-'^tree -^tree 

where Al^^^ and S^ree come from the tree level, and and B'^l^ come from the one- 
loop radiative corrections. The relationship between the hadronic radiative corrections of 
Eqs. (2.5-2.7) and the electron- quark radiative corrections of Eq. (2.16) follows from the 
linear combinations: 

Rl=' = {glRl-gt.RUl{9l-9v). 
RP = {9vRv + 9vRv + 9vRv)/i9v + 9v + 9v), (2-17) 

- -^A + ^Ay 

Ra~^ = (Ra + R'aiI'^i 

R^A ~ ~Ra "I" Ra "I" Ra- 

For the axial corrections we have used the explicit values of g\, which are +1 or —1, and 

the R^^^ correction is calibrated to 1, since the tree level term gA~^ 9a — ^■ 

Each particle carries 4-momentum p? = m? with the following structure in the center-of- 
mass (CM) frame: 

pi = (Ei,0,0,p), 

P2=(^2,0,0,-p), 

P3 = (E3,psin^,0,pcos^) , (2.18) 
P4 = (£^4, —p sin 9, 0, —p cos 9) . 

Here Ei = Eo^ and E2 = E^ (for elastic scattering) and p denote the energy and momentum 
of the scattered particles. The Mandelstam variables are defined as follows: 

s^{p^+P2f^{E^ + E2f^El^, 

t = (Pi -Psf = -2p' (1 - cos^) , (2.19) 
u^{pi- Pa? = (^1 - E2) - 2p2 (1 + cos 9) , 

with El 2 = p^ + m\2- For elastic scattering we have s-\-t-\-u — 2mi + 2m2. Input kinematic 
parameters, in our case, are the energy of the electron in the laboratory reference frame, 
E'lab, and the scattering angle 9. Hence Eqm and p can be expressed in terms of E'lab and 9: 



(Ecu + ^2 -^^If 
2E? 



P'- ' oz.2^ --mi (2.20) 



For quarks bound in the nucleon, we can assume that the energy of the quark can be defined 
as a fraction of the mass of the nucleon -Bfiab — x^'^'n^ with x ^ 1/3. Taking into account 
lab — P2,iab + "^2) ^c Can determine center-of-mass energy as follows: 

-^CM = (Pl.lab +P2,lab)^ 

/ {El,^^-miy/\{xmNy-miy/^ COS 9 A , ^ 

^ml + ml + 2E,^ub x 1 - ^ . 2.21 

V Ei^i^^xniN J 
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In the latter equation, 6i denotes the unknown angle between the spatial momentum of 
the electron and quark just before scattering took place. Final expressions for radiative 
corrections have been integrated over this angle to account for all possible values of 6i. 



III. DETAILS OF THE CALCULATION 
A. One-loop corrections 

Doing one-loop calculations in the Standard Model by hand is a laborious task due to the 
sheer number of particles and the complexity of the underlying theory. To facilitate this task, 
and to reduce the possibility of errors, for some time software packages have been developed 
to automate different tasks in these calculations [14]. In this paper, we have chosen to use 
the packages FeynArts, FormCalc, and LoopTools [9]. These packages are designed to work 
hand in hand on various aspects of the calculation, as described below. 

In total, 438 one-loop diagrams were calculated, including counterterms and contribu- 
tions from scalar bosons {H, x, 4>) and ghost fields {uz, For analysis and comparison 
with earlier work (see Ref. [4], for example), it is useful to split diagrams into three classes 
(Fig. 2). Class 1 are self-energy loops. Class 2 are triangles, and Class 3 are boxes. Class 
1 have a dominant contribution from Z exchange and jZ mixing. Class 2 include triangles 
with f f'Z and //'7 vertices. Class 3 are the ZZ, WW and 7Z box and crossed box di- 
agrams. Contributions from scalar bosons and ghost fields to the boxes were found to be 
negligibly small. 




FIG. 2. Classes of one-loop contributions to electron-quark scattering include self-energy, ver- 
tex, and box diagrams. 
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Because of the huge size of the analytical expressions for amplitudes, we leave them out 
of this article. Certain steps involved in the evaluation of one-loop PNC amplitude have to 
be explained more carefully. To evaluate four dimensional one-loop tensor integrals, we have 
used the Constrained Differential Renormalization (CDR) scheme [15] implemented in the 
package FormCalc. In the CDR scheme, Feynman diagrams are considered completely in four 
dimensions. Thereafter, the reduction of singular basic functions (products of propagators 
and their derivatives) has been renormalized into the sum of "regular" ones by implementing 
a set of rules in such a way that Ward identities are satisfied. It was proven in Ref. [9] that 
CDR is equivalent at the one- loop level to regularization by dimensional reduction [16], after 
taking the Fourier transform of the basic renormalized functions into momentum space. This 
last approach corresponds to a modified dimensional regularization, where one-loop integrals 
are considered in D dimensions, but all the tensors and spinors are kept 4-dimensional [16]. 
To preserve gauge invariance in dimensional reduction, one should use ^^j^ (jjj^ = with 
gn^g^'' = for the tensor decomposition. 

In dimensional regularization, the general structure of one-loop tensor integral can be 
written in the form [17] 

jN ^ (27r/x)^"^ f^D^ qnqi2---% ^3^^ 

Here, (yu) means the regularization scale parameter of dimensional reduction, which is related 
to the CDR renormalization scale by log ( ) = log {ji^) + 2. For the tensor decomposition 



of Eq. (3.1) into linear combinations of tensor coefficients functions, we have used the method 
of Passarino-Veltman [18]. 

The Passarino-Veltman approach deals with two, three, and four point tensor integrals 
with two, three, and four propagators, respectively. The rank of the tensors is equal to the 
number of integrable momenta [qi^qi2 • • • Q'jp^in the numerator of Eq. (3.1). Two, three and 
four point tensor coefficient functions, coming from tensor decomposition, were calculated 
numerically by using the package LoopTools, and algebraic calculations were completed with 
the help of packages FeynArts, FormCalc and FORM. To summarize this section, we provide 
an outline of how calculations were implemented in these packages. 



Generation of the diagrams (package FeynArts) 

• Topologies are defined 

• Fields are inserted 

• Amplitudes and their counterterms in integral form were defined, along with CKM 
matrix 

Evaluation of amplitudes (package FormCalc along with FORM) 

• Indices contracted and traces are taken 

• Amplitudes in general form (i.e. combination of many-point tensor coefficients along 
with spinor chains) were presented 
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• Radiative corrections were calculated (On Shell Renormalization) 

• Infrared divergences are treated by adding soft-photon emission 
Numerical evaluation 

• Standard Model parameters and kinematics are defined 

• Many-point tensor coefficients are numerically evaluated (package LoopTools) 

B. Renormalization constants 

Generally, tensor coefficient functions are ultraviolet divergent (inversely proportional to 

the parameter e = A — D). In order to cancel divergences and transform bare parameters 
into physical observables one has to introduce a renormalization scheme. The renormalized 
parameters are related to the bare parameters (denoted by a subscript 0) as follows: 

Co = (1 + 5e) e, 
(ZA _(l + \5Z^^ \5Z^^ \ ( Z\ 
\A^) -\ \bZ^^ 1 + \bZ^^ J I ^ J ' 

± 



(3.2) 



= (l + W 
Ho^[l + i^Z^) H, 

Counterterms were chosen in the On Shell Renormalization (OSR) scheme in the 't Hooft- 
Feynman gauge, where the gauge parameter ^ = 1. The renormalization constants are [19]: 

Wave function renormalization: 



5Z^^ = -Re ^Ef^ (e) , 5Z^^ = 2Re 

5Z^^ = -Re f ^Sl^ {k')] , = -2Re ^^^^ ^^^^ 



SZ-^-^ = -Re(^^Er{k^)) , SZ-^-Re(^^E-{k^)) , (3.3) 
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5Z^ = -Re 



'A. 



SZ'^ = -Re 



'_d_ 



SZli"^ = -Re (S^'^ (m| 



szl^^^ 



Re {ml 
_d_ 



-m^.Re 



Mass renormalization: 



5Ml = Re i E^^ ( M 



z ' 



5M^ = Re(s^^ (m^)), 



(3.4) 



= Re (e^ (m^)) , 5m/, = |m/,Re (m)) + E^'^ (mjj + 2E^'^ (m^ ) . 

Here L and 7? correspond to left- and right-handed fermions, E means one-loop integral of 
the truncated self-energy graph, and ± denotes the transverse part only. 
Charge and mixing angle renormalization: 



5 siv? 9w = cos^ 9 



5Ml 5M\ 



w 



w 



Ml 



5 cos^ Ow = —5 sin^ 9 



W: 



5e^--(sZ^^+'A^SZ'^' 
2 \ cos^^^y , 



(3.5) 



with sin^ = 1 - M^/M^. 



C. Infrared divergences and soft-photon emission 

For diagrams with photon exchange between external fermion legs we have encountered 
infrared (IR) divergences. These are regulated by introducing a small photon "rest mass" 
A in loop integrals involving the photon propagator. This results in In (m'j / A^^ terms in 
the amplitude. Hence the radiative corrections are, strictly speaking, infinite in the limit 
A — > 0. This unphysical dependence on A is cancelled by adding the inelastic bremsstrahlung 
contributions arising from soft-photon emission by the external fermions (i.e. from the initial 
or final electron or quark) to the scattering cross section [20]. 

To sec how this works in the parity-violating asymmetry, consider the inelastic 
bremsstrahlung amplitude Ai^, which includes the four diagrams of Fig. 3 involving Z ex- 
change, and the corresponding QED amplitude A41 (not shown) involving photon exchange. 
The total inelastic cross section is given by 



Ml +Mi-. (3.6) 
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FIG. 3. Bremsstrahlung diagrams for treating IR divergences. 



In analogy with Eq. (2.1), a parity- violating asymmetry will result from the interference 
term 2Re {{Ml)* M^} . 

In the soft photon approximation, it is straightforward to show (see Ref. [21], for example) 
that the amplitude for soft photon emission can be written in terms of the Born (tree level) 
amplitude 



Ml 



Born 



Pi 



P3 



QqP2 QqPA 



(pi-k) {p3-k) {p2-k) ipi-k) 



(3.7) 



with k the momentum of the photon and e{k) its polarization. This is just the amplitude 
for elastic scattering (without bremsstrahlung) times a factor for soft photon emission. An 
identical factor relates Ai^ to A^forn- 

The contribution of the parity- violating inelastic differential cross section to the radiative 
correction can also be expressed in terms of the parity-violating elastic Born cross section, 
once the appropriate integrals over phase space are done. The steps involved in carrying out 
the phase space integrals are analogous to ordinary parity-conserving electron scattering, as 
discussed for example in the paper by Maximon and Tjon [21]. We refer the reader to that 
paper for details, and simply present the result that 



where k^^^ is the soft-photon factor defined as 

(fk ( Pi , P3 



K. 



eq 
soft 



a 



/ 



+ 



+ 



QqP2 QqPi 



27r2 J\k\<Ae 2yFTA2 V (Pi • k) (p3 • k) (p2 • k) {p4 ■ k) 



(3.8) 



(3.9) 
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Here Ae is the maximum momentum of the undetected photon for which an elastic scattering 
event is recorded. It is related to the final electron detector acceptance in the lab frame, 

Adding together the elastic and inelastic cross sections, we have 



+ 



■MBorn + -Miorn 



eq 

soft' 



(3.10) 



where Ai"' 



and M.fg^ denote the one-loop radiative corrections. Note that at this stage, 
our expression (3.10) also includes the purely QED radiative corrections Al^ad to the Born 
amplitude A^Born- -^rad includes vertex and 77 box diagrams of order 0{a^). 

For the PV contribution, it is convenient to work with the amplitude Ai^y defined in 
Eq. (2.15), from which we can write the relevant PV interference term as 



da 



eq 

PV,tot 



jeq 

'PV,tree 



PV,rad 



+ (A^pV.tree)* (^^Jad + A^LrnlCft) } 



-'^'PV,tree2 



2 "-soft ) 



(3.11) 



At this point we note that the second term of Eq. (3.11) involves a PV interference term 
between A^^ad -Atpv.tree) pl^^^ the soft photon emission contribution needed to cancel the 
IR divergences in M^^- In principal, one should use the full expression (3.11) to evaluate 
both the QED and weak radiative corrections on the same footing, since they contribute 
to the PV asymmetry at the same order. In practice, this has not been done. We take 
the position here that experimental analyses have already accounted for the QED radiative 
corrections, and to include such effects again here would be double counting. Accordingly, 
we drop the second term in Eq. (3.11) from further consideration at this time. 

Using the general structure of the PV amplitude in Eq. (2.15), we can rewrite the first 
term in Eq. (3.11) as 



da^V,tot ~ 2Re (A1^_) 



eq 
VA 



_l_ _l_ P"? _L 1 L-'^l \ T^^l T 
"T l^-"- "T "T 2'^soft ) ^tree-J . 



eq 
AV 



(3.12) 



where current products Jyj^ and J^y are defined in Eq. (2.15). 
Evaluating Eq. (3.9), we have 



eq 
soft 



27r2 



2mll (pi,pi) - (2mi - t) I (^1,^3) + 2gJmJ/ (^2,^2) 
-Ql (2^g (P2, Pa) - 2Qq (u-ml- mfj I (pi, P4) 

-2Qq (s - - mfj I (pi,P2) 



(3.13) 



with the soft- photon emission integral I {pi,pj) defined as 

d^k 1 



\k\<Ae 2\/PTA2 {pi ■ k) {pj ■ k) 



(3.14) 
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These integrals generally have been worked out in Ref. [22]. For electron-quark scattering 
in the CM frame, we have (see also Ref. [19]) 



27raij 



2 2 2 



+ - In' 



'Ei-Pi 



Ei+pi 



--In' 



+U, l-^{E,+p,)]+uJl-^{E,- p,) 



' Ej-pA 
Ej+Pj) 

(3.15) 



-Lis 



^ {E^ + pM - Li2 il-^ {E, - pj) 



where Vij = {oif^ml — m^)/(2 {aijEi — Ej j), and {Ei,pi) are the energy and momentum in 
the CM system. The parameters for different values oi i,j are given in Table I. Li2 is the 
dilogarithm function 



^ [ dt 

J z 



ln(l - t) 



(3.16) 



TABLE I. Soft photon emission integral parameters of Eq. (3.15). 



r 

2 

1 

2 

1 
1 



1 

2 
3 
4 

4 

2 



rUe 
rriq 

rUe 

rriq 

rUe 



rriq 
rUe 
rriq 

rriq 
m„ 



1 

1 + 



2m? 



2mi 

9 



2m^ 



■mi —m^ ) — 4m^m^ 



■,2_ 



2mi 



s—mt—m: 



m|+m^— s) —Am'^m'^ 



The IR divergent (In A) terms in the elastic and inelastic contributions to the radiative 
correction cancel exactly. In particular, if we define the modified radiative correction factors 

Ry = Ry + I'^soft; (3.17a) 
R\ ^ R\ + \i%liv (3.17b) 

then Ry and R\ are the radiative corrections to electron-quark scattering free of IR diver- 
gences. This was confirmed numerically in our calculation by varying the parameter A over 
several orders of magnitude, with no numerically significant change in R. 

There remains a weak logarithmic dependence on As of the form ln(A£/m). This depen- 
dence will be cancelled if hard photon bremsstrahlung is taken into account. This requires 
knowledge of particular experimental details such as detector geometry, energy resolution. 
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phase space cuts, etc. This is beyond the scope of this paper, and we have chosen to simply 
set the scale for these effects by choosing an energy resolution AE" corresponding to the 
parameters of the SAMPLE detector. We leave a more complete treatment of the hard 
bremsstrahlung corrections, as well as the simultaneous treatment of both QED and weak 
radiative corrections, as a future project. 



IV. NUMERICAL RESULTS 

The parameters of the Standard Model are taken from Ref. [23], and arc given in Table II. 
As discussed by Hollik [24], the masses of the light quarks arc regarded as parameters, and 
are adjusted to reproduce the results of a dispersion analysis of the experimentally measured 
hadronic vacuum polarization. 

TABLE IL Standard model parameters used in this calculation. 



Quantity Value Quantity Value 

m'u 47. MeV rUe 0.51100 MeV 

m-d 47. MeV 105.66 MeV 

rUs 150. MeV 1777.0 MeV 

rUc 1.25 GeV Mz 91.1882 GeV 

nib 4.2 GeV Mw 80.419 GeV 

mt 174.3 GeV Mh 100. GeV 



Wc have chosen the SAMPLE kinematics, corresponding to ii^e,iab = 194 MeV and an 
average backward angle 9 — 146.2°. Radiative correction results for the V{e) x A{q) and 
A{e) X V{q) current interactions are given separately in the Table III for u, d, and s quarks. 

The contributions from self-energy, triangle, and box classes of diagrams are shown for 
loops which do not involve photons. The sclf-cncrgy (SE) diagrams give the dominant con- 
tribution to radiative corrections in most of the cases. Triangle and box diagrams that 
involve photon loops are shown as a combined result in column 5 because they are sepa- 
rately infrared divergent. In the combined result, the infrared divergence is cancelled by the 
bremsstrahlung corrections for soft-photon emission. The corrections show a weak logarith- 
mic dependence on the resolution of the detector, AE, which we have taken to be 120 MeV 
in this table. This corresponds roughly to the parameters of the SAMPLE detector, which 
detects all electrons above the Cerenkov threshold of 20 MeV [25] . 

The calculated one-quark radiative corrections were combined to form hadronic vector 
and axial vector corrections. These arc shown in Table IV for the isoscalar and isovector 
representations i?y^°, R\r^ ^ ^ ^i^d R^^- Also shown are the appropriate linear combi- 
nations of isoscalar and isovector corrections for the proton [Ry and R\) and the neutron 
[Ry and i?^). In addition, radiative corrections for the new SAMPLE kinematics, with 
£;iab = 120 MeV and denoted SAMPLE II, are also shown. Differences with the SAMPLE I 
results are mainly due to the bremsstrahlung corrections. 



15 



TABLE III. Radiative corrections for one-quark currents at SAMPLE I kinematics, 
(£;iab = 194 MeV, e = 146.2°, = 120 MeV). 



o//' 


Self Energy Vertex (no 7) 


Box (no 7) 


Vertex-|-Box (with 7) 


Total 


ue 
de 
se 
eu 
ed 
es 


-0.059 
-0.025 
-0.025 
-0.278 
-0.278 
-0.293 


-0.015 
0.026 
0.026 
-0.224 
-0.164 
-0.198 


0.027 
0.002 
0.002 
0.096 
0.022 
0.022 


0.019 
0.004 
0.015 
0.114 
-0.083 
-0.081 


-0.028 
0.008 
0.017 
-0.292 
-U.oUo 
-0.550 


TABLE IV. Radiative 
(Slab = 194 MeV, e = 146.2° 
e = 146.2°, = 77 MeV). 


corrections for nucleon currents at SAMPLE I 
, /S.E = 120 MeV), and SAMPLE II kinematics {Ei^b 


kinematics, 
= 120 MeV, 






T=0 


T=l 


isosinglet p 


n 


SAMPLE I Rv 

Ra 


0.057 
-0.210 


-0.005 

-0.398 


0.029 -0.252 
-0.760 -0.608 


0.022 
-0.187 


SAMPLE II Rv 

Ra 


0.053 
-0.216 


-0.012 
-0.408 


0.024 -0.271 
-0.781 -0.624 


0.017 
-0.192 



We point out again that these values are calculated in the CDR scheme, and also include 
a bremsstrahlung contribution. Values calculated in the MS scheme tend to be smaller [7] , 
however the tree level amplitudes also differ between the two schemes due to the different 
definitions of sin^ Oy^. 

It is interesting to examine the dependencies of the total radiative corrections -R^^^ and 
i?y on various kinematic and Standard Model parameters. These two corrections are large 
because the tree level amplitudes are suppressed by (1 — As\v?9w)- We have found only a 
weak dependence of the radiative corrections on scattering angle 6 and electron energy i?iab- 

Figure 4 shows the strong dependence of radiative corrections on the mass of the top 
quark. It also facilitates comparison with previous work [5], which used rrit = 120 MeV. 
This strong dependence on rrit arises from the large mass splitting within the top-bottom 
fermion doublet. The corrections show only a shght dependence on the Higgs mass Mh, 
which we have not shown. In Fig. 5 we show the (logarithmic) dependence of Ry and RjT^ 
on the the photon detection parameter AE. 

V. EXPERIMENT 

The first experiment to measure the weak neutral magnetic form factor of the nucleon 
is the SAMPLE experiment. The details of the techniques employed in this experiment are 
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(GeV) 

top ^ ' 

FIG. 4. -R^^^ (dashed line) and Ry (solid line) as a function of the mass of the top quark (in 
GeV). 



available in Refs. [10,12,25]. The experiment involves scattering 194 MeV polarized electrons 
from protons and deuterium. The scattered electrons are detected at 4 backward angles 
138° <Q < 160°, which results in an average ~ 0.1(GeV/c)2. 

The measured proton and deuteron asymmetries for Q"^ — 0.1 (GeV/c)^ and ^avg — 146.2° 
are (in ppm): 



^exp ^ _4_g2 ^ 0.61 ± 0.73, 
^^^P = -6.97 ±0.64 ±0.55. 



(5.1a) 
(5.1b) 



The first uncertainty is statistical while the second is systematic. In a subsequent analysis, 
Spayde [26] has reexamined the data, and finds two significant corrections. The first arises 
from a different treatment of the electromagnetic radiative corrections, while the second is 
a correction for the fraction of the signal due to pions. Each correction increases and 
by about 4 percent, resulting in quoted asymmetries of 



A;^p = -5.61 ±0.67 ±0.88, 
^exp ^ _7 28 ± 0.68 ± 0.75. 



(5.2a) 
(5.2b) 



Our theoretical results for the one quark contributions given in Table IV need to be 
supplemented by the many quark contributions. These have been evaluated using chiral 

perturbation theory by Zhu et al. [8], who found small, predominantly isovector, contribu- 
tions. Their isovector contribution can be parameterized as R^^^ = —0.06 ±0.24, while the 
isoscalar is Rj,^^ = —0.01 ± 0.14. The quoted uncertainties are large, and dominate the 
overall uncertainties in the theoretical analysis, as described below. 
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20 40 60 80 100 120 140 160 

AE (MeV) 

FIG. 5. i?5=^ (dashed hne) and (soUd hne) as a function of the photon detection parameter 
AE (in MeV). 

Our error analysis for the theoretical asymmetry takes into account uncertainties in the 
measured electromagnetic form factors, in the axial contributions G\ and G\, and in the 
radiative corrections. The largest source of uncertainty in the one quark radiative corrections 
is in the soft photon approximation for the bremsstrahlung contributions. We have allowed 
for a rather generous error estimate by evaluating these corrections for AE — 60 MeV, which 
is half the value used in Table IV. 

Our theoretical asymmetry for the SAMPLE kinematics, including all radiative correc- 
tions, and adding errors in quadrature, can be written to show the explicit dependence on 
the strange form factors: 

A^^ = (-7.29 ± 0.65) + (3.61 ± 0.07)^1^(0.1) + 1.94^^(0.1), (5.3a) 
= (-8.74 ± 0.89) + (0.82 ± 0.05)^^(0.1) + 1.46G'|;(0.1). (5.3b) 

For no strange quark contribution, the measured asymmetries are smaller than the theo- 
retical ones. The discrepancy between the theoretical and experimental values for Ap and 
Ad are about the same, whereas the coefficients of the dominant G%[ term are considerably 
different. 

In our analysis we have used the recent nucleon electromagnetic form factors determined 
by Brash et al. [11]. At = 0.1 GeV, these nucleon form factors tend to be smaller than 
in the Galster parameterization, thereby exacerbating the discrepancy between theoretical 
and experimental values noted in previous work [25,26]. The axial form factors are taken to 
be a simple dipole form, following Refs. [7,25]. 

The contribution from for these kinematic is expected to be small. For the range of 
values of the strange radius = ±0.22 found in the literature, we expect G|;(0.1) to be in the 
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range =p0.10. For comparison, the neutron electric form factor has the value (7^(0.1) = 0.036. 
HAPPEX [27] has measured the asymmetry Ap at 9 = 12.3° and = 0.477. They find the 
linear combination {G% + 0.392G%j) = 0.025 ± 0.020 ± 0.014, which would seem to favor a 
positive value of Gg. A measurement at = 0.1 is underway, and could shed some light 
on the SAMPLE result. 

Because the largest overall uncertainty in the theoretical calculations is from R^^^, it is 
useful to follow the analysis of the SAMPLE papers [25] and isolate this term explicitly. In 
this case, we have 

A^^ = (-5.93 ± 0.27) + (3.61 ± 0.07)G^(0.1) + (2.34 ± 0.03)G^'^=^(0.1), (5.4a) 
Af = (-7.10 ± 0.40) + (0.82 ± 0.05)^^(0.1) + (2.83 ± 0.04)G'^'^=^(0.1). (5.4b) 

We have dropped the G% terms, and have added (in quadrature) an additional uncertainty 
±0.18 to Ap and ±0.14 to A^. The isoscalar terms in G^, shown in Eq. (2.7) have been 
absorbed into the constant term. The axial isovector contribution is 

GY=\0.1) = -{1 + Rl=')Gl=\OA), (5.5) 

which has the theoretically expected value G^'"^~^(0.1) = —0.58 ± 0.27. This includes the 
large effect of the one quark axial correction (-R^^^ = —0.40 ±0.02), and the highly uncertain 
many quark contribution {Rl=^ = -0.06 ± 0.24) of Ref. [8]. 

Because the coefficients of G^' ~ in Eqs. (5.4) are roughly equal, changes in G^ ~ will 
shift Ap and A^ by about the same amount. This is consistent with the uniform discrepancy 
between theory and experiment in Ap and Ad referred to above, and so it is useful to try 
to fit the data in this way. However, there may be other explanations for such a uniform 
discrepancy, including other physics, or systematic errors in the experimental analysis of the 
kind determined by Spayde [26] . Hence we regard this as simply a device for characterizing 
the data, and indicating possible directions for explanation. 

Solving the pair of equations (5.4) for G%i and G^'^~^, and using the measured asym- 
metries from Eq. (5.2), we find 

G^(Q2 = 0.1) = 0.16 ±0.50, (5.6a) 
GY^^iQ'^ = 0.1) = -0.11 ± 0.49. (5.6b) 

Combining the errors in the extracted value of G^ and the theoretical value quoted above, 
the difference is about 1 a. 

There is some uncertainty about how to extrapolate G%j{Q'^ — 0.1) to jig — G%j{Q'^ = 0). 
A simple dipole form consistent with the proton and neutron magnetic electromagnetic form 
factors would give jis = 0.21 ± 0.65. Using the extrapolation model proposed by Hemmert 
et al. [28], one instead finds a value of the strange magnetic moment /x^ = 0.04 ±0.50 ±0.07, 
Here the last value of the error was introduced in the extrapolation. A negative value of 
is consistent with expectations from a number of theoretical models [29-32]. In particular, 
recent lattice QCD calculations find = -0.16 ± 0.18 [31] and /x, = -0.28 ± 0.10 [32]. By 
contrast, a recent SU(3) chiral quark-soliton model gives = 0.074 to 0.115 [33]. 

We have repeated the above analysis for the SAMPLE II kinematics {E = 120 MeV). 
This experiment is currently being analyzed. We find 
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A*^ = (-2.13 ± 0.08) + (1.07 ± 0.02)^^(0.043) + (1.06 ± 0.01)G^'^-^ (0.043), (5.7a) 
^ (_2 62 ± 0.13) + (0.26 ± 0.02)^^(0.043) + (1.35 ± 0.01)G'5'^=^(0.043), (5.7b) 

with the theoretically expected value (0.043) = —0.62 ± 0.28. The overall asymme- 

tries are smaller, and the uncertainties enter somewhat differently than in the SAMPLE I 
kinematics. This should help in determining the role of G^' ~ in the asymmetries. 

VI. SUMMARY 

We have reexamined the one-loop radiative corrections to the one-quark amplitudes in 
the CDR scheme. Our calculations include an evaluation of brcmsstrahlung corrections in 
the soft photon approximation. The brcmsstrahlung corrections cancel infrared divergences 
in the one-loop diagrams. The soft photon approximation has a logarithmic dependence on 
detection threshold. The sensitivity to this parameter will be reduced if hard bremsstrahlung 
photons are included, but that requires knowledge of the particular experimental setup. 

Although we find that G^' ~ is highly suppressed, there is no indication from the present 
calculation that it is suppressed by the magnitude indicated by the SAMPLE analysis. 
A contribution from the isosinglet strange axial form factor G% is unlikely to resolve the 
problem, as we find it is highly suppressed by the isosinglet radiative corrections i?^^ = 
—0.760. The many quark corrections calculated by Zhu et al. [8] were numerically small, 
but given the large uncertainties, this is one avenue for further theoretical work. Further 
experimental measurements by the HAPPEX collaboration at = 0.1 GcV^, as well as a 
new measurement by the SAMPLE group at the smaller electron energy 120 MeV, should 
also shed fight on this situation. 
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